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(57) A method for measuring the diameter of a micro 
liquid drop and the spatial distribution of micro liquid 
drops by measuring the diameter of an unsharp image 
due to defocusing and the number of interference fring- 
es thereof is applied to measurement of micro bubbles 
and applicable even if the spatial distribution density of 
micro liquid drops or micro bubbles is high. A liquid 
space where a micro bubble (10) is floating is irradiated 



with a sheet-like parallel laser beam (2), and a defocus 
image of the micro bubble (1 0) to which the laser beam 
is applied is picked up on a defocus plane (8) through 
an objective lens (6) in a direction at an angle 6 with 
respect to the direction where the laser beam travels. 
The interference fringes (9) in the defocus image (10") 
of the micro bubble (10) are counted, and the diameter 
of the micro bubble (10) is determined according to for- 
mula (4). 
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Description 

Technical Field 

[0001] The present invention relates to a method and 
apparatus, together with an optical system, for measur- 
ing the diameter, distribution and so forth of micro liquid 
droplets and micro gas bubbles. More particularly, the 
present invention relates to a method and apparatus, 
together with an optical system, for simultaneously 
measuring the diameter and distribution of micro liquid 
droplets and micro gas bubbles distributed in a space 
by an interferometric method. 

Background Art 

[0002] A method of accurately measuring the distribu- 
tion and diameters of micro liquid droplets of fuel inject- 
ed into an engine, for example, is demanded. Similarly, 
a method of accurately measuring the distribution and 
diameters of micro liquid droplets sprayed in the air is 
demanded to design a nozzle used in the spray dry 
method, for example. Further, a method of accurately 
measuring the diameter and distribution of gas bubbles, 
together with changes thereof, is demanded in the study 
of absorption of C0 2 in air bubbles into the sea and the 
behavior of gas bubbles in beer and wine. 
[0003] Thus, there is a strong demand in various fields 
for a method and apparatus for accurately measuring 
the diameter and distribution of micro liquid droplets and 
gas bubbles in the state of being present in a space. 
[0004] Regarding micro liquid droplets, there has 
heretofore been a method in which micro liquid droplets 
distributed in a space are photographed and the photo- 
graph is analyzed. This method involves a problem in 
terms of measurement accuracy because the photo- 
graph may be out of focus or may become unsharp for 
other reasons . The method further suffers from the prob- 
lem that real-time processing cannot be performed. A 
method in which the photograph is taken with a CCD 
camera is also known. This method also suffers from 
the problem in terms of measurement accuracy and the 
problem that real-time processing cannot be performed. 
Further, the method involves the problem that a great 
deal of time is required for analysis. A holographic tech- 
nique and a method using a CCD camera for imaging 
are also known. However, these methods similarly in- 
volve the problem in terms of measurement accuracy 
and the problems that real-time processing cannot be 
performed and a great deal of time is required for anal- 
ysis. There is also known a method in which the shad- 
ows of micro liquid droplets are captured directly with a 
CCD camera in order to obtain real-time capability. With 
this method, however, it is difficult to measure small par- 
ticles because of the influence of diffraction. The method 
further involves the problem that it is difficult to measure 
the diameter of micro liquid droplets in limited three-di- 
mensional positions. 



[0005] In addition, there has heretofore been known 
a method in which a plurality of particles are simultane- 
ously measured by specifying positions in a three-di- 
mensional space with a method known as LDV, PDA, 

s PDPA, etc. With this method, two laser beams are 
crossed in the air to form spatial interference fringes, 
and light scattered from liquid droplets crossing the in- 
terference fringes is observed with the same measure- 
ment volume from a plurality of different points. The di- 

'0 ameters of the micro liquid droplets are measured from 
the phase differences between the measurement sig- 
nals In this case, because the diameter of each individ- 
ual particle passing through the interference fringe area 
is measured, the method suffers from the problem that 

'5 measurement in the space surrounding the interference 
fringe area cannot simultaneously be performed. The 
measurement accuracy is also unsatisfactory. 
[0006] Under these circumstances, a method has 
been proposed (SAE Paper no. 950457, 960830) in 

20 which a sheet-shaped parallel laser beam is applied to 
a measurement space, and out-of-focus images of mi- 
cro liquid droplets irradiated with the laser beam are 
captured. In the out-of-focus image corresponding to 
each micro liquid droplet, interference fringes are 

is present, and there is a fixed relationship between the 
number of interference fringes present in the out-of-fo- 
cus image and the diameter of the micro liquid droplet. 
Accordingly, the diameter of the micro liquid droplet can 
be measured by measuring the number of interference 

30 fringes. It is also possible to measure the spatial distri- 
bution of the micro liquid droplets. 
[0007] With the above-described method of measur- 
ing the diameter and spatial distribution of micro liquid 
droplets by measuring the number of interference fring- 

35 es in each out-of-focus image, the applicable field is lim- 
ited to micro liquid droplets. The method has not here- 
tofore been applied to micro gas bubbles. 
[0008] Further, the above-described method involves 
the problem that when the spatial distribution density of 

40 micro liquid droplets is high, out-of-focus images over- 
lap each other because they are circular and occupy 
large areas. Therefore, it is difficult to measure the di- 
ameters of the micro liquid droplets separately. 



[0009] The present invention was made in view of the 
above-described problems with the prior art, and an ob- 
ject of the present invention is to expand the method of 
measuring the diameter and spatial distribution of micro 
liquid droplets by measuring the diameter of each out- 
of-focus image obtained by defocusing and the number 
of interference fringes in the out-of-focus image into a 
method of measuring the diameter and spatial distribu- 
tion of micro gas bubbles, and to provide a measuring 
optical system that allows the method to be applied to 
a case where the spatial distribution density of micro liq- 
uid droplets and micro gas bubbles is high. 
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[0010] Another object of the present invention is to 
provide a method and apparatus for determining the po- 
sition, diameter and velocity of micro liquid droplets and 
micro gas bubbles from the analysis of out-of -focus im- 
ages. 

[0011] A method of measuring the diameter, distribu- 
tion and so forth of micro gas bubbles according to the 
present invention, which is provided to attain the above- 
described objects, is a method wherein a sheet-shaped 
parallel laser beam is applied to a liquid space in which 
micro gas bubbles are floating, and out-of-focus images 
of micro gas bubbles irradiated with the laser beam are 
captured from a lateral direction which is at an angle 9 
to the direction of travel of the laser beam. The number 
N of interference fringes in the out-of-focus image cor- 
responding to each micro gas bubble is measured, and 
the diameter D of the micro gas bubble is determined 
from the following relationship: 

D=(2X.N/na) [cos (8/2)-sin (6/2) 

+V{n 2 +1-2ncos(9/2)}]' 1 (4) 

where X is the wavelength of the laser beam; a is 
the angle subtended at the micro gas bubble by an ob- 
jective lens used to capture the image of the micro gas 
bubble; and n is the relative index of refraction of a liquid 
in which the micro gas bubble is present. 
[0012] Another method of measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 
liquid droplets according to the present invention is a 
method wherein a sheet-shaped parallel laser beam is 
applied to a space in which micro gas bubbles or micro 
liquid droplets are floating; out-of-focus images of micro 
gas bubbles or micro liquid droplets irradiated with the 
laser beam are captured from a lateral direction which 
is at a predetermined angle to the direction of travel of 
the laser beam; and the numbers of interference fringes 
in the respective out-of-focus images corresponding to 
the micro gas bubbles or the micro liquid droplets are 
measured to determine the diameters and distribution 
of the micro gas bubbles or the micro liquid droplets. 
[0013] The method is characterized in that the out-of- 
focus images are captured with an imaging optical sys- 
tem at an imaging plane where the images are out of 
focus in a direction parallel to a plane containing the di- 
rection of travel of the sheet-shaped parallel laser beam 
and an optical axis of the imaging optical system and 
where the images are substantially in focus in a direction 
perpendicular to the plane. 

[0014] In this case, it is desirable that the spacing of 
interference fringes on the imaging plane be adjustable 
by adjusting the defocus condition of the out-of-focus 
images. 

[0015] The arrangement may be such that the sheet- 
shaped parallel laser beam is moved in parallel to a di- 



4 

rection perpendicular to the plane of the sheet-shaped 
parallel laser beam with respect to the space in which 
micro gas bubbles or micro liquid droplets are floating, 
and the out-of-focus images are captured in synchro- 
5 nism with the movement of the sheet-shaped parallel 
laser beam 

[0016] Another method of measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 
liquid droplets according to the present invention is a 

10 method wherein a sheet-shaped parallel laser beam is 
applied to a space in which micro gas bubbles or micro 
liquid droplets are floating, and linear out-of-focus im- 
ages of micro gas bubbles or micro liquid droplets irra- 
diated with the laser beam are captured from a lateral 

'5 direction which is at a predetermined angle to the direc- 
tion of travel of the laser beam with an imaging optical 
system at an imaging plane where the images are out 
of focus in a direction parallel to a plane containing the 
direction of travel of the laser beam and an optical axis 

20 of the imaging optical system and where the images are 
substantially in focus in a direction perpendicular to the 
plane. The linear out-of-focus images extend in the di- 
rection of the plane in correspondence to the micro gas 
bubbles or the micro liquid droplets. The center of each 

25 of the out-of-focus images is determined to thereby de- 
termine the center position of the corresponding micro 
gas bubble or micro liquid droplet. 
[0017] In this case, it is desirable that the center po- 
sition be determined from a peak position of a moving 

30 average value obtained by taking an average in the 
range extending from a distance L/2 forward of a specific 
position to a distance L/2 rearward of the specific posi- 
tion in the longitudinal direction and determining the av- 
erage to be a value at this position, where L is the length 

35 of a linear out-of-focus image, and successively moving 
the specific position. 

[0018] Another method of measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 
liquid droplets according to the present invention is a 

40 method wherein a sheet-shaped parallel laser beam is 
applied to a space in which micro gas bubbles or micro 
liquid droplets are floating, and linear out-of-focus im- 
ages of micro gas bubbles or micro liquid droplets irra- 
diated with the laser beam are captured from a lateral 

is direction which is at a predetermined angle to the direc- 
tion of travel of the laser beam with an imaging optical 
system at an imaging plane where the images are out 
of focus in a direction parallel to a plane containing the 
direction of travel of the laser beam and an optical axis 

so of the imaging optical system and where the images are 
substantially in focus in a direction perpendicular to the 
plane. The linear out-of-focus images extend in the di- 
rection of the plane in correspondence to the micro gas 
bubbles or the micro liquid droplets. Each of the out-of- 

55 focus images is subjected to Fourier transform to obtain 
a frequency, and the obtained frequency is multiplied by 
the length of the out-of-focus image to obtain the 
number of interference fringes in the out-of-focus image. 
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The diameter of the micro gas bubble or the micro liquid 
droplet is determined on the basis of the number of in- 
terference fringes. 

[0019] In this case, it is desirable that discrete Fourier 
transform be performed as the Fourier transform to ob- 
tain a discrete frequency distribution, and function fitting 
be applied to the discrete frequency distribution to ob- 
tain the diameter of the micro gas bubble or the micro 
liquid droplet. 

[0020] A further method of measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 
liquid droplets according to the present invention is a 
method wherein a sheet-shaped parallel laser beam is 
applied to a space in which micro gas bubbles or micro 
liquid droplets are floating, and two image frames are 
captured at a micro time interval At, each of which two 
image frames contains linear out-of-focus images of mi- 
cro gas bubbles or micro liquid droplets irradiated with 
the laser beam. The linear out-of-focus images are cap- 
tured from a lateral direction which is at a predetermined 
angle to the direction of travel of the laser beam with an 
imaging optical system at an imaging plane where the 
images are out of focus in a direction parallel to a plane 
containing the direction of travel of the laser beam and 
an optical axis of the imaging optical system and where 
the images are substantially in focus in a direction per- 
pendicular to the plane. The linear out-of-focus images 
extend in the direction of the plane in correspondence 
to the micro gas bubbles or the micro liquid droplets. 
Cross correlation between the two captured image 
frames is calculated for each linear out-of-focus image 
in the two captured image frames to obtain the displace- 
ment As, of each linear out-of-focus image, and the ve- 
locity U| of each micro gas bubble or micro liquid droplet 
is determined from the following relationship: 

u,=As/At (6) 

[0021] In this case, it is desirable to remove a high- 
frequency component corresponding to interference 
fringes in the linear out-of-focus image when calculating 
cross correlation between the two captured image 
frames. 

[0022] A still further method of measuring the diame- 
ter, distribution and so forth of micro gas bubbles and 
micro liquid droplets according to the present invention 
is a method wherein a sheet-shaped parallel laser beam 
is applied to a space in which micro gas bubbles or micro 
liquid droplets are floating, and linear out-of-focus im- 
ages of micro gas bubbles or micro liquid droplets irra- 
diated with the laser beam are captured from a lateral 
direction which is at a predetermined angle to the direc- 
tion of travel of the laser beam with an imaging optical 
system at an imaging plane where the images are out 
of focus in a direction parallel to a plane containing the 
direction of travel of the laser beam and an optical axis 
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of the imaging optical system and where the images are 
substantially in focus in a direction perpendicular to the 
plane. The linear out-of-focus images extend in the di- 
rection of the plane in correspondence to the micro gas 
5 bubbles or the micro liquid droplets. The center of each 
of the out-of-focus images is determined to thereby de- 
termine the center position of the corresponding micro 
gas bubble or micro liquid droplet. Each of the out-of- 
focus images is subjected to Fourier transform to obtain 
10 afrequency, and the obtained frequency is multiplied by 
the length of the out-of-focus image to obtain the 
number of interference fringes in the out-of-focus image. 
The diameter of the micro gas bubble or the micro liquid 
droplet is determined on the basis of the number of in- 
'5 terference fringes. Further, two image frames containing 
the linear out-of-focus images are captured at a micro 
time interval At. Cross correlation between the two cap- 
tured image frames is calculated for each linear out-of- 
focus image in the two captured image frames to obtain 
20 the displacement As, of each linear out-of-focus image, 
and the velocity u, of each micro gas bubble or micro 
liquid droplet is determined from the following relation- 
ship: 

U|=As,/At (6) 

[0023] An apparatus for measuring the diameter, dis- 
tribution and so forth of micro gas bubbles and micro 

30 liquid droplets according to the present invention in- 
cludes laser beam application means for applying a 
sheet-shaped parallel laser beam to a space in which 
micro gas bubbles or micro liquid droplets are floating, 
and imaging means for capturing linear out-of-focus im- 

35 ages of micro gas bubbles or micro liquid droplets irra- 
diated with the laser beam, which is applied by the laser 
beam application means, from a lateral direction which 
is at a predetermined angle to the direction of travel of 
the laser beam with an imaging optical system at an im- 

40 aging plane where the images are out of focus in a di- 
rection parallel to a plane containing the direction of trav- 
el of the laser beam and an optical axis of the imaging 
optical system and where the images are substantially 
in focus in a direction perpendicular to the plane. The 

45 linear out-of-focus images extend in the direction of the 
plane in correspondence to the micro gas bubbles or the 
micro liquid droplets. The apparatus further includes 
center position measuring means for determining the 
center of each of the out-of-focus images to thereby de- 

50 termine the center position of the corresponding micro 
gas bubble or micro liquid droplet, and diameter meas- 
uring means for subjecting each of the out-of-focus im- 
ages to Fourier transform to obtain a frequency, multi- 
plying the obtained frequency by the length of the out- 

55 of-focus image to obtain the number of interference 
fringes in the out-of-focus image, and determining the 
diameter of the micro gas bubble or the micro liquid 
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droplet on the basis of the number of interference fring- 
es. Further, the apparatus includes velocity measuring 
means for capturing two image frames containing the 
linear out-of-focus images at a micro time interval At, 
calculating cross correlation between the two captured 
image frames for each linear out-of-focus image in the 
two captured image frames to obtain the displacement 
AS| of each linear out-of-focus image, and determining 
the velocity u, of each micro gas bubble or micro liquid 
droplet from the following relationship: 

UpAs/At (6) 

[0024] An optical system for measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 
liquid droplets according to the present invention is a 
measuring optical system wherein a sheet-shaped par- 
allel laser beam is applied to a space in which micro gas 
bubbles or micro liquid droplets are floating; out-of-fo- 
cus images of micro gas bubbles or micro liquid droplets 
irradiated with the laser beam are captured from a lateral 
direction which is at a predetermined angle to the direc- 
tion of travel of the laser beam; and the numbers of in- 
terference fringes in the respective out-of-focus images 
corresponding to the micro gas bubbles or the micro liq- 
uid droplets are measured to determine the diameters 
and distribution of the micro gas bubbles or the micro 
liquid droplets. 

[0025] The measuring optical system is characterized 
by including an imaging optical system in which the focal 
length or the image-side principal plane in a direction 
parallel to a plane containing the direction of travel of 
the sheet-shaped parallel laser beam and an optical axis 
of the imaging optical system and the focal length or the 
image-side principal plane in a direction perpendicular 
to the plane containing the optical axis of the imaging 
optical system are different from each other, and image 
pickup means placed in an image plane which is in the 
vicinity of the image-formation plane in the direction per- 
pendicular to the above-described plane and which is 
off the image-formation plane in the direction parallel to 
the above-described plane. 

[0026] In this case, it is desirable that the imaging op- 
tical system be an anamorphic optical system compris- 
ing a combination of an axially symmetric objective lens 
and a cylindrical lens. 

[0027] It is also desirable that at least one of the focal 
length and the image-side principal plane of the imaging 
optical system in the direction parallel to the plane be 
adjustable. 

[0028] It is also desirable that the imaging optical sys- 
tem have a rectangular aperture elongated in the direc- 
tion parallel to the plane. 

[0029] In the present invention, out-of-focus images 
of micro gas bubbles or micro liquid droplets are cap- 
tured with an imaging optical system at an imaging plane 



where the images are out of focus in a direction parallel 
to a plane containing the direction of travel of a sheet- 
shaped parallel laser beam and the optical axis of the 
imaging optical system and where the images are sub- 

5 stantially in focus in a direction perpendicular to the 
plane. Consequently, the out-of-focus image corre- 
sponding to each micro gas bubble or micro liquid drop- 
let becomes a one-dimensional image compressed in 
the direction perpendicular to the plane. Therefore, even 

10 when the spatial distribution density of micro gas bub- 
bles and micro liquid droplets is high, the respective out- 
of-focus images can be separated from each other. Ac- 
cordingly, the number of interference fringes in each out- 
of-focus image can be readily counted separately from 

'5 each other In addition, it becomes easy to determine 
the center position of each out-of-focus image to detect 
the distributed conditions of micro gas bubbles or micro 
liquid droplets. Even in such a case, the position, diam- 
eter and velocity distributions of micro gas bubbles and 

20 micro liquid droplets can be measured simultaneously 
and accurately. 

Brief Description of the Drawings 
25 [0030] 

Fig. 1 is a diagramfor explaining the principle ofthe 
method of measuring the diameter and spatial dis- 
tribution of micro gas bubbles according to the 
30 present invention and the principle of a convention- 
al method of measuring the diameter and spatial 
distribution of micro liquid droplets and shows an 
example of an out-of-focus image of a micro gas 
bubble or a micro liquid droplet. 
35 Fig. 2 is a diagram for analyzing light rays passing 
through a micro gas bubble floating in a liquid. 
Fig. 3 is a diagram for analyzing light rays passing 
through a micro liquid droplet floating in the air. 
Fig. 4 is a perspective view showing a first embod- 
40 iment of the optical system for measuring the diam- 
eter, distribution and so forth of micro gas bubbles 
and micro liquid droplets according to the present 
invention. 

Fig. 5 is a perspective view showing a second em- 
45 bodiment of the optical system for measuring the 
diameter, distribution and so forth of micro gas bub- 
bles and micro liquid droplets according to the 
present invention. 

Fig. 6 is a diagram showing an example of out-of- 
so focus images captured with an arrangement as 
shown in Fig. 1(a). 

Fig. 7 is a diagram showing out-of-focus images 
captured with an arrangement as shown in Fig. 4, 
the out-of-focus images corresponding to those 
55 shown in Fig. 6 

Fig. 8 is a diagram showing an example of an image 
signal obtained with regard to one interference 
fringe image. 
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Fig. 9 is a diagram showing image signals of inter- 
ference fringe images obtained at an image pickup 
surface and an example of the result of taking mov- 
ing average. 

Fig. 10 is a diagram showing an image signal of an 
interference fringe image and an example of a pow- 
er spectrum obtained therefrom by FFT. 
Fig. 1 1 is a diagram for explaining a method of ac- 
curately obtaining the frequency of the original sig- 
nal from a discrete power spectrum by function fit- 
ting. 

Fig. 12 is a diagram for explaining a method of de- 
termining the velocity of micro liquid droplets by cal- 
culating cross correlation. 
Fig. 13 is a diagram showing an example of an im- 
age frame displaying the position, diameter and ve- 
locity distributions of micro liquid droplets or micro 
gas bubbles measured simultaneously. 

Best Mode for Carrying Out the Invention 

[0031] The following is a description of the principles 
and embodiments of the method and apparatus for 
measuring the diameter, distribution and so forth of mi- 
cro liquid droplets and micro gas bubbles and the meas- 
uring optical system according to the present invention. 
[0032] First, the principle of a publicly known method 
of measuring the diameter and spatial distribution of mi- 
cro liquid droplets by measuring the number of interfer- 
ence fringes in each out-of-focus image will be de- 
scribed with a view to facilitating understanding. 
[0033] First, as shown in Fig. 3, a plane wave 2 is 
made incident on a micro liquid droplet 1 of refractive 
index n floating in the air. At this time, there is a differ- 
ence in the angle 6 measured when twice-refracted light 
4 of incident angle (hereinafter, both the incident angle 
and the refraction angle will be assumed to be angles 
measured from a tangential plane to the interface) i-, 
and once-reflected light 3 of incident angle % 0 are paral- 
lel to each other and the phase difference therebetween 
is 2rTm (m is an integer) and when the phase difference 
between the refracted light 4 and the reflected light 3 is 
2(m+1)7c. The angular difference A6 is given by 

A9=(2X/D)[n sin(6/2) 
+V {n 2 +1-2n cos(e/2)}+cos(e/2)]" 1 (1) 

where 6 is the viewing angle of scattered light from 
the micro liquid droplet 1 with respect to the illuminating 
light 2; D is the diameter of the micro liquid droplet 1; 
and X is the wavelength of the illuminating light 2. 
[0034] The meaning of the above is as follows. As 
shown in Fig. 1 (a), scattered light 5 from the micro liquid 
droplet 1 contains a row of high-intensity portions (inter- 
ference fringes) produced by interference at a micro an- 



gular spacing A0, which are centered in the direction of 
scattering angle 9 with respect to the illuminating light 
2. When an objective lens (imaging lens) 6 is placed in 
the path of the scattered light 5 to form an image V of 

5 the micro liquid droplet 1 on an image plane 7 by the 
scattered light 5, an out-of-focus image 1" of the micro 
liquid droplet 1 such as that shown in Fig 1(b) is ob- 
tained on an out-of-focus plane (defocus plane) 8, which 
is off the image plane 7. The range shown by the dashed 

10 lines in Figs. 1 (a) and (b) indicates the range of a bundle 
of rays incident on the objective lens 6. The external size 
and shape of the out-of-focus image 1 " of the micro liq- 
uid droplet 1 obtained on the out-of-focus plane 8 de- 
pend on the size of the objective lens 6 and the distance 

'5 from the image plane 7 to the out-of-focus plane 8 inde- 
pendently of the size of the micro liquid droplet 1 . When 
the external shape of the objective lens 6 is circular, the 
out-of-focus image 1 " of the micro liquid droplet 1 is cir- 
cular. The number N of interference fringes 9 formed 

20 within the circle is determined by the angle a subtended 
at the micro liquid droplet 1 by the objective lens 6 and 
the above-described angular difference A9. 
[0035] That is, from the relationship of o=NxA6 and 
the above equation (1), the diameter of the micro liquid 

25 droplet 1 is given by 

D= (2XN/a) [n sin(9/2) 
30 +V<n 2 +1-2ncos(9/2)}+cos(9/2)]" 1 (2) 

[0036] The diameter D of the micro liquid droplet 1 can 
be obtained by substituting the number N of interference 
fringes 9 in the out-of-focus image 1 " actually observed 

35 and measured into equation (2). 

[0037] As will be clear from Fig. 1 (a), when sheet- 
shaped parallel light extending in a direction perpendic- 
ular to the plane of the figure is used as the illuminating 
light 2 and micro liquid droplets 1 1 2 , ■■■ are present in 

40 the path of the light in addition to the micro liquid droplet 
1, out-of-focus images 1,", 1 2 ", - are also obtained on 
the out-of-focus plane 8 as in the case of the micro liquid 
droplet 1 , and the diameter D can similarly be obtained. 
The center positions of the out-of-focus images 1 ,", 1 2 ", 

45 - approximately correspond to the center positions of 
the images 1', 1,', 1 2 ', ••• of the micro liquid droplets 1 , 
1,, 1 2 , - on the image plane 7 Therefore, the distribu- 
tion of the micro liquid droplets and the diameter of each 
micro liquid droplet can be simultaneously determined 

50 from the out-of-focus images 1,", 1 2 ", ■•• obtained on the 
out-of-focus plane 8. 

[0038] The foregoing is the principle of the publicly 
known method of measuring the diameter and spatial 
distribution of micro liquid droplets by measuring the 
55 number of interference fringes in each out-of-focus im- 
age. Let us consider obtaining the distribution and di- 
ameter of micro gas bubbles present in a liquid in place 
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of micro liquid droplets. 
[0039] In Fig. 2, a plane wave 2 is made incident on 
a micro gas bubble 10 of refractive index 1 floating in a 
liquid. At this time, there is a difference in the angle 6 
measured when twice-refracted light 12 of incident an- 
gle t, and once-reflected light 1 1 of incident angle i 0 are 
parallel to each other and the phase difference therebe- 
tween is 2rrm (m is an integer) and when the phase dif- 
ference between the refracted light 1 2 and the reflected 
light 1 1 is 2(m+1 )7t. The angular difference A6 is given by 

A8=(2X/nD) [cos (8/2)-sin (6/2) 

+V{n 2 +1-2ncos(9/2)}]' 1 (3) 

where 8 is the viewing angle of scattered light from 
the micro gas bubble 1 0 with respect to the illuminating 
light 2; D is the diameter of the micro gas bubble 1 0; and 
X is the wavelength of the illuminating light 2. 
[0040] The meaning of the above is as follows. As 
shown in Fig. 1(a), scattered light 5 from the micro gas 
bubble 10 contains a row of high-intensity portions (in- 
terference fringes) produced by interference at a micro 
angular spacing A8, which are centered in the direction 
of scattering angle 6 with respect to the illuminating light 
2. When an objective lens 6 is placed in the path of the 
scattered light 5 to form an image 10' of the micro gas 
bubble 1 0 on an image plane 7 by the scattered light 5, 
an out-of-focus image 10" of the micro gas bubble 10 
such as that shown in Fig. 1 (b) is obtained on an out-of- 
focus plane (defocus plane) 8, which is off the image 
plane 7. The range shown by the dashed lines in Figs. 
1 (a) and (b) indicates the range of a bundle of rays in- 
cident on the objective lens 6. The external size and 
shape of the out-of-focus image 10" of the micro gas 
bubble 1 0 obtained on the out-of-focus plane 8 depend 
on the size of the objective lens 6 and the distance from 
the image plane 7 to the out-of-focus plane 8 independ- 
ently of the size of the micro gas bubble 10. When the 
external shape of the objective lens 6 is circular, the out- 
of-focus image 1 0" of the micro gas bubble 1 0 is circular. 
The number N of interference fringes 9 formed within 
the circle is determined by the angle a subtended at the 
micro gas bubble 10 by the objective lens 6 and the 
above-described angular difference AG. 
[0041] That is, from the relationship of a=NxA9 and 
the above equation (3), the diameter D of the micro gas 
bubble 10 is given by 

D=(2XN/noc) [cos (9/2)-sin(8/2) 

+V {n 2 +1-2ncos (8/2))]" 1 (4) 

[0042] The diameter D of the micro gas bubble 1 0 can 
be obtained by substituting the number N of interference 
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fringes 9 in the out-of-focus image 1 0" actually observed 
and measured into equation (4). 
[0043] As will be clear from Fig. 1(a), when sheet- 
shaped parallel light extending in a direction perpendic- 
5 ular to the plane of the figure is used as the illuminating 
light 2 and micro gas bubbles 10,, 10 2 , - are present 
in the path of the light in addition to the micro gas bubble 
10, out-of-focus images 10,", 10 2 ", - are also obtained 
on the out-of-focus plane 8 as in the case of the micro 
io gas bubble 1 0, and the diameter D can similarly be ob- 
tained The center positions of the out-of-focus images 
10,", 01 2 ", — approximately correspond to the center 
positions of the images 10'. 10,', 10 2 ', ■■• of the micro 
gas bubbles 10, 10,, 10 2 , •■• on the image plane 7. 
'5 Therefore, the distribution of the micro gas bubbles and 
the diameter of each micro gas bubble can be simulta- 
neously determined from the out-of-focus images 10,", 
10 2 ", •■■ obtained on the out-of-focus plane 8. 
[0044] It will be understood from the above discussion 
20 that in the case of micro gas bubbles also, when a sheet- 
shaped parallel laser beam is applied to a measurement 
space to capture out-of-focus images of micro gas bub- 
bles irradiated with the laser beam, interference fringes 
are present in the out-of-focus image corresponding to 
25 each micro gas bubble, and there is a fixed proportional 
relationship between the number of interference fringes 
present in the out-of-focus image and the diameter of 
the micro gas bubble. Accordingly, it is possible to meas- 
ure the diameter of the micro gas bubble by measuring 
30 the number of interference fringes, and the distribution 
of the micro gas bubbles can be simultaneously ob- 
tained from the distribution of the center positions of the 
out-of-focus images 

[0045] Incidentally, Fig. 6 shows schematically an ex- 
35 ample of out-of-focus images captured with an arrange- 
ment as shown in Fig. 1 (a) when the spatial distribution 
density of micro liquid droplets or micro gas bubbles is 
high In the following, we will discuss micro liquid drop- 
lets typically because it is understood that micro gas 
40 bubbles and micro liquid droplets can be handled simi- 
larly except for the difference between equations (4) and 
(2). 

[0046] Fig. 6 shows out-of-focus images a, b, c and d 
of four micro liquid droplets 1 captured at the out-of-fo- 
« cus plane 8 with the arrangement shown in Fig. 1 (a) in 
a case where the four micro liquid droplets 1 are present 
in close proximity to each other in the path of the sheet- 
shaped parallel illuminating light 2. Because the four mi- 
cro liquid droplets 1 are too close to each other, the out- 
50 of-focus images a, b, c and d with circular outer shapes 
corresponding to the micro liquid droplets 1 overlap 
each other. Accordingly, it is not easy to count the 
number of interference fringes 9 in each of the images 
a, b, c and d separately from each other. It is also difficult 
55 to detect the distributed conditions of micro liquid drop- 
lets 1 by determining the center positions of the images 

[0047] Accordingly, an optical system as shown in the 
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perspective view of Fig. 4 is used as a first embodiment 
of the optical system for measuring the diameter, distri- 
bution and so forth of micro gas bubbles and micro liquid 
droplets according to the present invention. First, let us 
define a coordinate system. The direction of travel of 
sheet-shaped parallel illuminating light 2 to be applied 
to micro liquid droplets 1 , 1 1 , 1 2 , ••• is denoted by S, and 
the optical axis of a measuring optical system 20 is de- 
noted by O. The optical axis 0 is set in a plane perpen- 
dicular to the plane of the sheet-shaped parallel light 2. 
The direction perpendicular to the optical axis O in that 
plane is defined as an x-axis direction, and the direction 
perpendicular to both the optical axis 0 and the x-axis 
direction and parallel to the sheet-shaped parallel illu- 
minating light 2 is defined as a y-axis direction. The 
measuring optical system 20 shown in Fig. 4 includes 
an objective lens 6 and a cylindrical lens 21 (a negative 
cylindrical lens in the case of Fig. 4) placed in coaxial 
relation to the objective lens 6 and having a refracting 
power only in the x-axis direction (having no refracting 
power in the y-axis direction). An image pickup surface 
22 of an image pickup device, e.g. a CCD, is placed in 
the image-formation plane in the y-axis direction of the 
measuring optical system 20, that is, in the image-for- 
mation plane of the objective lens 6. In contrast, the im- 
age-formation plane in the x-axis direction of the meas- 
uring optical system 20 is formed at a position off the 
image pickup surface 22 (behind the image pickup sur- 
face 22 in the case of Fig. 4), With this arrangement, the 
out-of-focus image of the micro liquid droplet 1 , which 
is located in the vicinity of the optical axis 0, for example, 
is circular in the optical path from the objective lens 6, 
which has a circular aperture, to the cylindrical lens 21 . 
However, as the distance from the cylindrical lens 21 
increases toward the image pickup surface 22, the out- 
of-focus image gradually increases in flatness, and on 
the image pickup surface 22, the out-of-focus image is 
a horizontal line. However, there is no change in the 
number of interference fringes 9 in the out-of-focus im- 
age at any position. 

[0048] Fig. 7 shows out-of-focus images a, b, c and d 
of four micro liquid droplets 1 and so forth obtained from 
the image pickup surface 22 with the arrangement 
shown in Fig. 4, which correspond to those in Fig. 6. It 
should be noted, however, that the each out-of-focus im- 
age is illustrated on the assumption that there is no 
change in the magnifications in the x- and y-axis direc- 
tions (in reality, because the focal length in the x-axis 
direction and so forth may change, the magnification of 
the out-of-focus image may also change). As will be 
clear from the comparison of Figs. 6 and 7, the out-of- 
focus images a, b, c and d captured with the arrange- 
ment shown in Fig. 4 are those which are obtained by 
compressing the circular out-of-focus images a, b, c and 
d captured with the arrangement shown in Fig. 1(a) in 
the vertical direction (y-axis direction) with the center 
positions thereof left as they are, thereby converting 
them into one-dimensional images (in the x-axis direc- 



tion). Accordingly, the four out-of-focus images a, b, c 
and d no longer overlap each other in the y-axis direc- 
tion, so that the number of interference fringes 9 in each 
of the images a, b, c and d can be counted separately 

5 with ease. It also becomes easy to determi ne the center 
positions of the images a, b, c and d to thereby detect 
the distributed conditions of the micro liquid droplets 1 
and so forth (this will be described later). 
[0049] It should be noted that the out-of-focus images 

to a, b, c and d as shown in Fig. 6, which are captured by 
using an axially symmetric measuring optical system, 
have circular edges around them. Therefore, the diam- 
eter of each image can easily be found and it is easy to 
count the number of interference fringes 9 in the aper- 

'5 ture. In the case of the out-of-focus images a, b, c and 
d compressed as in Fig. 7, however, the light quantity at 
the center of each image is large. Consequently, the 
light quantity in the vicinity of each end is relatively 
small. Accordingly, the ends of the out-of-focus image 

20 are inconspicuous, and the length L thereof is unclear. 
However, if the measuring optical system is under the 
same conditions and the out-of-focus plane is the same, 
the length L of any compressed out-of-focus image re- 
mains the same. Therefore, there will be no problem in 

25 this regard if confirmation is made once in advance un- 
der the same conditions. 

[0050] It should be noted that compressing an out-of- 
focus image in the vertical direction (y-axis direction) al- 
lows an improvement in contrast of the captured out-of- 

30 focus image and enables the measurement sensitivity 
to increase, advantageously. 
[0051] Incidentally, the lens arrangement of the meas- 
u ring optical system 20 that produces an in-focus con- 
dition in the y-axis direction and an out-of-focus condi- 

35 tion in the x-axis direction on the image pickup surface 
22 as shown in Fig. 4 may be an anamorphic optical 
system comprising a combination of an axially symmet- 
ric objective lens 6 and a cylindrical lens 21 as stated 
above, and it may also be an anamorphic optical system 

40 using a plane-symmetry anamorphic surface, e.g. atoric 
surface, as a refracting surface. It is also possible to use 
an optical system in which the refracting power in the x- 
axis direction and that in the y-axis direction are the 
same but the principal plane in the x-axis direction and 

45 that in the y-axis direction are different from each other 
and which is therefore focused on the image pickup sur- 
face 22 in the y-axis direction but defocused at the im- 
age pickup surface 22 in the x-axis direction. The above- 
described optical systems may be arranged to include 

so a reflecting surface, as a matter of course. 

[0052] Fig. 5 is a perspective view of an optical system 
for measuring the diameter, distribution and so forth of 
micro gas bubbles and micro liquid droplets according 
to a second embodiment of the present invention, which 

55 is arranged to further ameliorate unsatisfactory points 
of the optical system shown in Fig. 4. The measuring 
optical system 20' comprises an objective lens 6 and a 
combination of a positive cylindrical lens 21 , and a neg- 
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ative cylindrical lens 21 2 which are placed in coaxial re- 
lation to the objective lens 6 and have a refracting power 
only in the x-axis direction (having no refracting power 
in the y-axis direction). The position of each of the two 
cylindrical lenses 21, and 21 2 is adjustable along the 
optical axis O. An image pickup surface 22 of an image 
pickup device is placed in the image-formation plane of 
the objective lens 6, which is the image-formation plane 
in the y-axis direction of the measuring optical system 
20'. 

[0053] With the above-described arrangement, the 
image-formation plane in the x-axis direction of the en- 
tire measuring optical system 20' can be adjusted freely 
with respect to the image pickup surface 22 by adjusting 
the relative positions of the positive cylindrical lens 21 , 
and the negative cylindrical lens 21 2 and the position of 
the combination of these cylindrical lenses relative to 
the objective lens 6. The focal length of the measuring 
optical system 20' in the x-axis direction can also be con- 
tinuously adjusted freely within a certain range. Accord- 
ingly, out-of-focus images (Fig. 7) compressed in the 
vertical direction (y-axis direction) into one-dimensional 
images (x-axis direction) are captured at the image pick- 
up surface 22 as in the case of Fig. 4. In addition, the 
length L of each linear out-of-focus image extending in 
the x-axis direction can be adjusted by adjusting the po- 
sitions of the two cylindrical lenses 21, and 21 2 
[0054] In the case of Fig. 7, the problem of overlap 
between the out-of-focus images a, b, c and d in the y- 
axis direction in Fig. 6 is resolved. However, out-of-focus 
images located at the same height (the same y-coordi- 
nate value) may overlap each other at their edge por- 
tions. Therefore, the overlap in the x-axis direction can- 
not be eliminated with the arrangement shown in Fig. 4. 
In such a case, if the arrangement shown in Fig. 5 is 
used, the overlap between the edge portions can be 
eliminated by making an adjustment so that the length 
L of each out-of-focus image is shortened. In this case 
also, there is no change in the number of interference 
fringes 9 in one out-of-focus image, as has been stated 
above. 

[0055] Further, as will be clear from the form of equa- 
tions (2) and (4), there is a proportional relationship be- 
tween the number N of interference fringes and the di- 
ameter D of a micro liquid droplet (micro gas bubble). 
Therefore, when the diameter D of each micro liquid 
droplet 1 under measurement is large, the number of 
interference fringes 9 in one out-of-focus image is large. 
Consequently, the interference fringes 9 in the captured 
image frame may become so fine that it is not easy to 
count the number of interference fringes 9. In such a 
case, the positions of the two cylindrical lenses 21 , and 
21 2 are adjusted so as to increase the length L of each 
out-of-focus image in reverse relation to the above, 
thereby increasing the resolution and thus making it 
possible to facilitate the counting of interference fringes 
9. 

[0056] Incidentally, in the arrangement shown in Fig. 



5, a slit-shaped aperture 23 elongated in the x-axis di- 
rection is placed in the vicinity of the objective lens 6 to 
limit the numerical aperture in the y-axis-direction so as 
to increase the depth of focus (depth of field). As a result, 
5 even when the optical axis O of the measuring optical 
system 20' is at an angle other than 90° with respect to 
the sheet-shaped parallel illuminating light 2, it is possi- 
ble to capture and measure an out-of-focus image of a 
micro liquid droplet 1 , , for example, which is somewhat 
'0 away from the optical axis O. It should be noted that be- 
cause the slit-shaped aperture 23 has a shape elongat- 
ed in the x-axis direction as stated above, it has no in- 
fluence on the number of interference fringes at a micro 
angular spacing AO that can be taken in for measure- 
's ment. Thus, the slit-shaped aperture 23 exerts no influ- 
ence on the number N of interference fringes in each 
individual out-of-focus image captured. 
[0057] Incidentally, as has been suggested above, the 
angle e of the optical axis O of the measuring optical 
20 system 20 or 20' with respect to the sheet-shaped par- 
allel illuminating light 2 is normally set at an angle be- 
tween 0° and 90". In this case, if the principal plane of 
the objective lens 6 and the image pickup surface 22 
are set at right angles to the optical axis O, it is difficult 
25 to image all micro liquid droplets in the oblique object 
plane 2 in the desired state unless the above-described 
slit-shaped aperture 23 is used. Therefore, a swing & tilt 
technique is adopted to tilt the principal plane of the ob- 
jective lens 6 and the image pickup surface 22 with re- 
30 spect to the optical axis 0 or to move them vertically by 
combining together shift, tilt and swing as used in pho- 
tography, thereby allowing all micro liquid droplets in the 
oblique object plane 2 to be imaged in the desired state. 
Examples of this technique include a method in which 
35 the principal plane of the objective lens 6 and the image 
pickup surface 22 are tilted with respect to the optical 
axis O so as to satisfy shine proof conditions. 
[0058] In the foregoing, the sheet-shaped parallel il- 
luminating light 2 is applied to the measurement space 
40 to determine the distribution and diameter of micro liquid 
droplets or micro gas bubbles located in the illumination 
sheet plane. However, it is possible to determine the dis- 
tribution and diameter of micro liquid droplets or micro 
gas bubbles in a three-dimensional space by moving the 
« sheet-shaped parallel illuminating light 2 in a direction 
perpendicular to the illumination sheet plane and cap- 
turing out-of-focus images at the image pickup surface 
22 individually in synchronism with the movement of the 
sheet-shaped parallel illuminating light 2. In this case, it 
50 is preferable to move the image pickup surface 22 in the 
direction of the optical axis in association with the move- 
ment of the sheet-shaped parallel illuminating light 2. 
[0059] Let us further describe embodiments of a 
method and apparatus for determining the position, di- 
ss ameter and velocity of micro liquid droplets or micro gas 
bubbles by using out-of-focus images as shown in Fig. 
7, which are captured at the image pickup surface 22 of 
the above-described measuring optical system 20 ac- 
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cording to the present invention. 
[0060] From the above discussion, the center posi- 
tions of the linear out-of-focus images a, b, c and d of 
length L obtained from the image pickup surface 22 ap- 
proximately correspond to the center positions of the mi- 5 
cro liquid droplets 1 -, , 1 2 , ••• in the sheet-shaped parallel 
light 2. Therefore, a method of determining the centers 
of the out-of-focus images (hereinafter referred to as "in- 
terference fringe images") a, b, c and d corresponding 
to the micro liquid droplets will be described first. A cap- '0 
tured image frame A as shown in Fig. 7 is horizontally 
scanned along the lengthwise direction (x-direction) of 
linear interference fringe images a, b, c and d and ver- 
tically scanned in a direction (y-direction) perpendicular 
to the lengthwise direction, thereby obtaining an image >5 
signal for the whole captured image frame A. The image 
signal contains a signal corresponding to each of the 
interference fringe images a, b, c and d as shown in Fig, 
8 by way of example. The signal has a length L in terms 
of distance and has peaks corresponding to the number 20 
N of interference fringes therein. Such signals are 
present in correspondence to the respective positions 
of the interference fringe images a, b, c and d. In Fig. 8, 
the abscissa axis corresponds to the position in the lon- 
gitudinal direction of the interference fringe image (ex- 25 
pressed in pixels), and the ordinate axis corresponds to 
the signal intensity. To determine the center position of 
a single interference fringe image as shown in Fig. 8, 
moving average should be taken along the longitudinal 
direction of the interference fringe image. The length L 30 
of one interference fringe image is determined by the 
condition of the measuring optical system 20 and the 
position of the image pickup surface 22. Therefore, the 
average is taken in the range extending from a distance 
L/2 forward of a specific position to a distance L/2 rear- 35 
ward of the specific position and determined to be a val- 
ue at this position, and the specific position is moved 
successively to obtain a moving average value. Fig. 9 
shows the image signals of the interference fringe im- 
ages obtained at the image pickup surface 22 and an <*o 
example of the result of taking moving average as stated 
above. When moving average is taken, an approximate- 
ly triangular-wave signal is obtained as illustrated in the 
figure. The peak position (1) of the signal is the center 
position of the interference fringe image corresponding *s 
to a micro liquid droplet. It should be noted that the range 
(3) in Fig. 9 corresponds to the length L of the interfer- 
ence fringe image 

[0061] Incidentally, a signal longer than the length L 
of the signal of an interference fringe image may appear 50 
in the signal in the same horizontal scanning direction. 
This occurs in some rare cases where a plurality of in- 
terference fringe images superimposed on one another 
are present in the same horizontal scanning direction. 
In this case, the half-width of the image signal of the 55 
interference fringe image or the half-width of the moving 
average signal becomes longer than in normal cases. 
Therefore, the situation can be readily judged. In such 



a case, no problem arises even if the image signal of 
the interference fringe image is removed. It is also pos- 
sible to determine the centers of two interference fringe 
images from the half -width. 

[0062] It should be noted that because a noise may 
be mixed in the image signal, it is desirable to judge that 
an interference fringe image is present only when it is 
decided that the amplitude (2) of the high-frequency 
component in Fig. 9 is more than a predetermined value. 
[0063] Thus, the position of each micro liquid droplet 
in the captured image frame A is determined, and the 
micro liquid droplet distribution and density in the space 
are determined. 

[0064] Next, a method of determining the diameter of 
each micro liquid droplet will be described. As has been 
stated above, because the length of a signal indicating 
each interference fringe image in the image signal ob- 
tained by scanning the captured image frame A is L, the 
range extending from L/2 forward of the center position 
determined as stated above to L/2 rearward of the cent- 
er position corresponds to the signal of each interfer- 
ence fringe image. Therefore, a signal extending over 
the range of L/2 forward of the determined center of the 
interference fringe image to L/2 rearward of it, that is, a 
signal having a length L in the longitudinal direction and 
centered at the center position of the interference fringe 
image, is cut out, and the absolute value or square (pow- 
er spectrum) of Fourier transform of the cut-out signal 
is determined, thereby obtaining the frequency f of the 
interference fringe image. By multiplying the frequency 
f by the length L of the interference fringe image, the 
number N of interference fringes in the out-of-focus im- 
age is obtained. Thereafter, N is substituted into equa- 
tion (2) or (4) to obtain the diameter D of the micro liquid 
droplet or the micro gas bubble. 
[0065] To Fourier-transform the signal of length L and 
determine the frequency f from the power spectrum , ac- 
tually, the signal is multiplied by a Hanning window func- 
tion, for example, as a window function for eliminating 
the influence of the edge, thereby performing fast Fou- 
rier transform (FFT). Incidentally, FFT is a kind of dis- 
crete Fourier transform. In discrete Fourier transform, 
the frequency interval obtained is 1/MA, where A is the 
sampling interval for a signal to be transformed and M 
is the sampling number. Accordingly, frequency can be 
obtained only at discrete frequency positions. If the fre- 
quency of the interference fringe image is precisely co- 
incident with any one of the frequencies discrete at a 
frequency interval 1/MA, the Fourier-transformed fre- 
quency signal appears as a single peak at the frequency 
position. However, when the frequency of the interfer- 
ence fringe image is present between two adjacent dis- 
crete frequencies, the signal undesirably appears not 
only at the positions of the two adjacent frequencies but 
also at the positions of discrete frequencies around 
them. An example of this is shown in Fig. 1 0. Fig. 1 0(a) 
shows a signal of an interference fringe image. Fig. 10 
(b) shows a power spectrum obtained by multiplying the 
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signal by a Harming window function and performing 
FFT. As will be clear from Fig. 10(b), a peak P K is present 
at the position of a frequency k, and signals P k .-, and 
P k+1 are also present at the positions of discrete fre- 
quencies k-1 and k+1 on both sides of the frequency k. 5 
There are also signals on both sides of these signals. 
To accurately obtain the frequency of the original signal 
by function fitting from such a discrete power spectrum, 
methods using various fitting functions are available. 
However, the frequency of the original signal can be ac- '0 
curately obtained by a method using a Gaussian func- 
tion (R.J. Adrian, et al. "Applications of Laser Tech- 
niques to Fluid Mechanics 5th International Symposium 
Lisbon, Portugal, 9-12 July, 1990" pp. 268-287 (Spring- 
er-Verlag)). That is, as shown in Fig. 11, when a peak <5 
P k is present at the position of a discrete frequency k 
and signals P k-1 and P k+1 smaller than the peak P k are 
also present at the positions of discrete frequencies k- 
1 and k+1 on both sides of the peak P k , the frequency f 
of the original signal can be obtained without consider- 20 
ing signals at other frequencies as follows: 



[0066] As the above-described window function and 
fitting function, other publicly known functions are also 
usable, as a matter of course. 30 
[0067] Next, a method of determining the velocity 
(vector) of each micro liquid droplet will be described. In 
this case, image frames A and A' as shown in Fig. 7 are 
captured at a micro time interval At. Assuming that the 
two captured image frames A and A' are those which 35 
are schematically shown in Figs. 1 2(a) and (b), interfer- 
ence fringe images a, b, c and d in each of the captured 
image frames A and A' change as illustrated in the fig- 
ures. Therefore, the cross correlation between the inter- 
ference fringe images a, b, c and d in the two capture 40 
image frames A and A 1 is calculated to determine the 
displacement As, of each interference fringe image in 
the form of a vector. The displacement distribution is 
schematically shown in Fig. 12 (c). 
[0068] From the displacement As, thus determined, 45 
the velocity u, of each micro liquid droplet is determined 
as follows: 

UpAs/At (6) 50 

[0069] More specifically, sheet-shaped parallel illumi- 
nating light 2 (Fig. 4 or 5) is applied at a micro time in- 
terval At by using a double-pulsed laser, for example, 
and two images A and A' are captured at the image pick- 55 
up surface 22 synchronously with the emission of the 
light. Signals of interference fringe images are cut out 



from the images A and A'. The method of cutting out 
signals is the same as the above-described method of 
determining the frequency. In this case, however, not 
only a signal from a single scanning line but also signals 
from adjacent scanning lines perpendicularly intersect- 
ing the interference fringe image are cut out at the same 
time, and the cross correlation is calculated for each in- 
terference fringe image cut out from the images A and 
A'. In the calculation of the cross correlation, each inter- 
ference fringe image cut out from the two images is shift- 
ed by one pixel at a time in the x- and y-directions to 
calculate a correlation value. The upper limit of the dis- 
placement for correlation calculation in the x- and y-di- 
rections is set appropriately in advance. The displace- 
ment of the vector to a position (peak position) where 
the highest correlation value is obtained is determined 
to be displacement As, for each interference fringe im- 
age. 

[0070] In the above-described calculation, each cut- 
out interference fringe image includes the interference 
fringe signal, and the interference fringes may move left- 
ward or rightward in the interference fringe image ac- 
cording to the phase. Therefore, the displacement of the 
interference fringe image obtained from the above-de- 
scribed cross correlation is not always the same as the 
actual displacement of the interference fringe image. 
Accordingly, if the cross correlation is calculated with re- 
spect to the interference fringe image in the state of in- 
cluding the interference fringe signal as stated above, it 
is not always possible to obtain an accurate displace- 
ment of the interference fringe image. Therefore, it is 
desirable to remove the high-frequency component 
from the interference fringe signal by passing the inter- 
ference fringe image signal through a low-pass filter be- 
fore the cross correlation is calculated. 
[0071] Further, because the calculation of the cross 
correlation is also performed in units of one pixel, an 
displacement As, in subpixel units cannot directly be ob- 
tained. A peak position expressed in subpixel units can 
be accurately obtained by applying various functions to 
discrete cross-correlation values obtained in the x- and 
y-directlons (e.g. a Gaussian function or a quadratic 
function may be used; however, a sine function or a co- 
sine function is preferably used because the cross cor- 
relation with a sine function is calculated in the present 
invention). 

[0072] It should be noted, however, that the displace- 
ment As, obtained by the above-described technique 
does not always accurately correspond to the displace- 
ment of a micro liquid droplet between different instanc- 
es in time. Therefore, it is desirable to judge whether or 
not the displacement As, accurately correspond as fol- 
lows. The frequency in the cut-out signal of each inter- 
ference fringe image is determined as stated above, and 
it is judged whether or not the frequency has changed. 
Alternatively, it is judged whether or not the rate of 
change due to evaporation or condensation is less than 
a predetermined value. Only when the rate of change is 



f=f k +1/2x{(logP k . r logP k+1 ) 

+ (logP k . r 2logP k +logP k+1 )} (5) 



(6) 5D 
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less than the predetermined value, it is judged that the 
displacement As, accurately corresponds to the dis- 
placement of the micro liquid droplet. 
[0073] Fig. 13 shows an example of an image frame 
displaying the position, diameter and velocity distribu- 5 
tions of micro liquid droplets or micro gas bubbles simul- 
taneously measured as stated above. The center of 
each circle indicates position. The size of each circle 
indicates diameter. The line segment indicates velocity. 
[0074] It should be noted that an apparatus for carry- >o 
ing out the above-described method of determining the 
position, diameter and velocity of micro liquid droplets 
or micro gas bubbles using out-of-focus images can be 
readily implemented via software using a personal com- 

[0075] Although the method and apparatus for meas- 
uring the diameter, distribution and so forth of micro gas 
bubbles and micro liquid droplets, together with the op- 
tical system for measuring the diameter, distribution and 
so forth of micro gas bubbles and micro liquid droplets, 20 
according to the present invention have been described 
above on the basis of embodiments, the present inven- 
tion is not limited to these embodiments but can be mod- 
ified in a variety of ways, 

25 

Industrial Applicability 

[0076] As will be clear from the foregoing description, 
the method and apparatus for measuring the diameter, 
distribution and so forth of micro gas bubbles and micro 30 
liquid droplets, together with the optical system for 
measuring the diameter, distribution and so forth of mi- 
cro gas bubbles and micro liquid droplets, according to 
the present invention are arranged to capture out-of-fo- 
cus images of micro gas bubbles or micro liquid droplets 
with an imaging optical system at an imaging plane 
where the images are out of focus in a direction parallel 
to a plane containing the direction of travel of a sheet- 
shaped parallel laser beam and the optical axis of the 
imaging optical system and where the images are sub- 
stantially in focus in a direction perpendicular to the 
plane. Consequently, the out-of-focus image corre- 
sponding to each micro gas bubble or micro liquid drop- 
let becomes a one-dimensional image compressed in 
the direction perpendicular to the plane. Therefore, even 
when the spatial distribution density of micro gas bub- 
bles and micro liquid droplets is high, the respective out- 
of-focus images can be separated from each other. Ac- 
cordingly, the number of interference fringes in each out- 
of-focus image can be readily counted separately from 
each other. In addition, it becomes easy to determine 
the center position of each out-of-focus image to detect 
the distributed conditions of micro gas bubbles or micro 
liquid droplets. Even in such a case, the position, diam- 
eter and velocity distributions of micro gas bubbles and 
micro liquid droplets can be measured simultaneously 
and accurately. 



A method of measuring a diameter, distribution and 
so forth of micro gas bubbles, said method compris- 
ing the steps of: 

applying a sheet-shaped parallel laser beam to 
a liquid space in which micro gas bubbles are 
floating; 

capturing out-of-focus images of micro gas 
bubbles irradiated with the laser beam from a 
lateral direction which is at an angle 9 to a di- 
rection of travel of the laser beam; 
measuring the number N of interference fringes 
in the out-of-focus image corresponding to 
each micro gas bubble; and 
determining a diameter D of the micro gas bub- 
ble from the following relationship: 

D = (2XN/na) [cos (6/2)-sin (6/2) 

H-V{n 2 +1-2ncos(e/2)}]" 1 (4) 

where X is a wavelength of the laser beam; a 
is an angle subtended at the micro gas bubble by 
an objective lens used to capture the image of the 
micro gas bubble; and n is a relative index of refrac- 
tion of a liquid in which the micro gas bubble is 
present. 

2. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets, said method comprising the steps of: 

applying a sheet-shaped parallel laser beam to 
a space in which micro gas bubbles or micro 
liquid droplets are floating; 
capturing out-of-focus images of micro gas 
bubbles or micro liquid droplets irradiated with 
the laser beam from a lateral direction which is 
at a predetermined angle to a direction of travel 
of the laser beam; and 

measuring numbers of interference fringes in 
the respective out-of-focus images corre- 
sponding to the micro gas bubbles or the micro 
liquid droplets to determine diameters and dis- 
tribution of the micro gas bubbles or the micro 
liquid droplets; 

wherein said out-of-focus images are cap- 
tured with an imaging optical system at an imaging 
plane where the images are out of focus in a direc- 
tion parallel to a plane containing the direction of 
travel of said sheet-shaped parallel laser beam and 
an optical axis of the imaging optical system and 
where the images are substantially in focus in a di- 
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rection perpendicular to said plane. 

3. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets according to claim 2, wherein a spacing of in- 
terference fringes on the imaging plane is adjusta- 
ble by adjusting a defocus condition of said out-of- 
focus images. 

4. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets according to claim 2 or 3, wherein said sheet- 
shaped parallel laser beam is moved in parallel to 
a direction perpendicular to a plane of said sheet- 
shaped parallel laser beam with respect to the 
space in which micro gas bubbles or micro liquid 
droplets are floating, and said out-of-focus images 
are captured in synchronism with movement of said 
sheet-shaped parallel laser beam. 

5. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets, said method comprising the steps of: 

applying a sheet-shaped parallel laser beam to 
a space in which micro gas bubbles or micro 
liquid droplets are floating; 
capturing linear out-of-focus images of micro 
gas bubbles or micro liquid droplets irradiated 
with the laser beam from a lateral direction 
which is at a predetermined angle to a direction 
of travel of the laser beam with an imaging op- 
tical system at an imaging plane where the im- 
ages are out of focus in a direction parallel to a 
plane containing the direction of travel of said 
laser beam and an optical axis of the imaging 
optical system and where the images are sub- 
stantially in focus in a direction perpendicular 
to said plane, said linear out-of-focus images 
extending in a direction of said plane in corre- 
spondence to the micro gas bubbles or the mi- 
cro liquid droplets; and 

determining a center of each of said out-of-fo- 
cus images, thereby determining a center po- 
sition of the corresponding micro gas bubble or 
micro liquid droplet. 

6. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets according to claim 5, wherein said center posi- 
tion is determined from a peak position of a moving 
average value obtained by taking an average in a 
range extending from a distance L72 forward of a 
specific position to a distance L/2 rearward of the 
specific position in a longitudinal direction and de- 
termining the average to be a value at this position, 
where L is a length of a linear out-of-focus image, 
and successively moving the specific position. 
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7. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets, said method comprising the steps of: 

5 applying a sheet-shaped parallel laser beam to 

a space in which micro gas bubbles or micro 
liquid droplets are floating; 
capturing linear out-of-focus images of micro 
gas bubbles or micro liquid droplets irradiated 

'0 with the laser beam from a lateral direction 

which is at a predetermined angle to a direction 
of travel of the laser beam with an imaging op- 
tical system at an imaging plane where the im- 
ages are out of focus in a direction parallel to a 

'5 plane containing the direction of travel of said 

laser beam and an optical axis of the imaging 
optical system and where the images are sub- 
stantially in focus in a direction perpendicular 
to said plane, said linear out-of-focus images 

20 extending in a direction of said plane in corre- 

spondence to the micro gas bubbles or the mi- 
cro liquid droplets; 

subjecting each of the out-of-focus images to 
Fourier transform to obtain a frequency; 
25 multiplying the obtained frequency by a length 

of the out-of-focus image to obtain the number 
of interference fringes in the out-of-focus im- 
age; and 

determining a diameter of the micro gas bubble 
so or the micro liquid droplet on a basis of the 

number of interference fringes. 

8. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 

35 lets according to claim 7, wherein discrete Fourier 
transform is performed as said Fourier transform to 
obtain a discrete frequency distribution, and func- 
tion fitting is applied to the discrete frequency dis- 
tribution to obtain a diameter of the micro gas bub- 

*o ble or the micro liquid droplet. 

9. A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets, said method comprising the steps of: 

45 

applying a sheet-shaped parallel laser beam to 
a space in which micro gas bubbles or micro 
liquid droplets are floating; 
capturing two image frames at a micro time in- 

50 terval At, the two image frames each containing 

linear out-of-focus images of micro gas bubbles 
or micro liquid droplets irradiated with the laser 
beam, the linear out-of-focus images being 
captured from a lateral direction which is at a 

55 predetermined angle to a direction of travel of 

the laser beam with an imaging optical system 
at an imaging plane where the images are out 
of focus in a direction parallel to a plane con- 
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taining the direction of travel of said laser beam 
and an optical axis of the imaging optical sys- 
tem and where the images are substantially in 
focus in a direction perpendicular to said plane, 
said linear out-of-focus images extending in a 
direction of said plane in correspondence to the 
micro gas bubbles or the micro liquid droplets; 
calculating cross correlation between the two 
captured image frames for each linear out-of- 
focus image in the two captured image frames 
to obtain a displacement As, of each linear out- 
of-focus image; and 

determining a velocity u { of each micro gas bub- 
ble or micro liquid droplet from the following re- 
lationship: 



of interference fringes in the out-of-focus im- 
age; 

determining a diameter of the micro gas bubble 
or the micro liquid droplet on a basis of the 
number of interference fringes; 
capturing two image frames containing the lin- 
ear out-of-focus images at a micro time interval 
At; 

calculating cross correlation between the two 
captured image frames for each linear out-of- 
focus image in the two captured image frames 
to obtain a displacement As, of each linear out- 
of-focus image; and 

determining a velocity u { of each micro gas bub- 
ble or micro liquid droplet from the following re- 
lationship: 



10. A method of measuring a diameter, distribution and 20 
so forth of micro gas bubbles and micro liquid drop- 
lets according to claim 9, wherein when cross cor- 
relation is calculated between said two captured im- 
age frames, a high-frequency component corre- 
sponding to interference fringes in the linear out-of- 25 
focus image is removed. 



1 1 . A method of measuring a diameter, distribution and 
so forth of micro gas bubbles and micro liquid drop- 
lets, said method comprising the steps of: 30 

applying a sheet-shaped parallel laser beam to 
a space in which micro gas bubbles or micro 
liquid droplets are floating; 
capturing linear out-of-focus images of micro 35 
gas bubbles or micro liquid droplets irradiated 
with the laser beam from a lateral direction 
which is at a predetermined angle to a direction 
of travel of the laser beam with an imaging op- 
tical system at an imaging plane where the im- 40 
ages are out of focus in a direction parallel to a 
plane containing the direction of travel of said 
laser beam and an optical axis of the imaging 
optical system and where the images are sub- 
stantially in focus in a direction perpendicular 45 
to said plane, said linear out-of-focus images 
extending in a direction of said plane in corre- 
spondence to the micro gas bubbles or the mi- 
cro liquid droplets; 

determining a center of each of said out-of-fo- 50 
cus images, thereby determining a center po- 
sition of the corresponding micro gas bubble or 
micro liquid droplet; 

subjecting each of the out-of-focus images to 
Fourier transform to obtain a frequency; 55 
multiplying the obtained frequency by a length 
of the out-of-focus image to obtain the number 



12. An apparatus for measuring a diameter, distribution 
and so forth of micro gas bubbles and micro liquid 
droplets, said apparatus comprising: 

laser beam application means for applying a 
sheet-shaped parallel laser beam to a space in 
which micro gas bubbles or micro liquid drop- 
lets are floating; 

imaging means for capturing linearout-of-focus 
images of micro gas bubbles or micro liquid 
droplets irradiated with the laser beam, which 
is applied by said laser beam application 
means, from a lateral direction which is at a pre- 
determined angle to a direction of travel of the 
laser beam with an imaging optical system at 
an imaging plane where the images are out of 
focus in a direction parallel to a plane contain- 
ing the direction of travel of said laser beam and 
an optical axis of the imaging optical system 
and where the images are substantially in focus 
in a direction perpendicular to said plane, said 
linear out-of-focus images extending in a direc- 
tion of said plane in correspondence to the mi- 
cro gas bubbles or the micro liquid droplets; 
center position measuring means for determin- 
ing a center of each of said out-of-focus images 
to thereby determine a center position of the 
corresponding micro gas bubble or micro liquid 
droplet; 

diameter measuring means for subjecting each 
of the out-of-focus images to Fourier transform 
to obtain a frequency, multiplying the obtained 
frequency by a length of the out-of-focus image 
to obtain the number of interference fringes in 
the out-of-focus image, and determining a di- 
ameter of the micro gas bubble or the micro liq- 
uid droplet on a basis of the number of interfer- 
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ence fringes; and 

velocity measuring means for capturing two im- 
age frames containing the linear out-of-focus 
images at a micro time interval At, calculating 
cross correlation between the two captured im- 
age frames for each linear out-of-focus image 
in the two captured image frames to obtain a 
displacement As, of each linear out-of-focus im- 
age, and determining a velocity u, of each micro 
gas bubble or micro liquid droplet from the fol- 
lowing relationship: 



at least one of the focal length and the image-side 
principal plane of said imaging optical system in the 
direction parallel to said plane is adjustable. 

16. An optical system for measuring a diameter, distri- 
bution and so forth of micro gas bubbles and micro 
liquid droplets according to any one of claims 12 to 
14, wherein said imaging optical system has a rec- 
tangular aperture elongated in the direction parallel 
to said plane. 



UpAS/At (6) 

15 

13. An optical system for measuring a diameter, distri- 
bution and so forth of micro gas bubbles and micro 
liquid droplets by applying a sheet-shaped parallel 
laser beam to a space in which micro gas bubbles 

or micro liquid droplets are floating, capturing out- 20 
of-focus images of micro gas bubbles or micro liquid 
droplets irradiated with the laser beam from a lateral 
direction which is at a predetermined angle to a di- 
rection of travel of the laser beam, and measuring 
numbers of interference fringes in the respective 25 
out-of-focus images corresponding to the micro gas 
bubbles or the micro liquid droplets to determine di- 
ameters and distribution of the micro gas bubbles 
or the micro liquid droplets, said optical system 
comprising: 30 

an imaging optical system in which a focal 
length or an image-side principal plane in a di- 
rection parallel to a plane containing the direc- 
tion of travel of said sheet-shaped parallel laser 35 
beam and an optical axis of the imaging optical 
system and a focal length or an image-side 
principal plane in a direction perpendicular to 
said plane containing the optical axis of the im- 
aging optical system are different from each 40 
other; and 

image pickup means placed in an image plane 
which is in a vicinity of an image-formation 
plane in the direction perpendicular to said 
plane and which is off an image-formation 45 
plane in the direction parallel to said plane. 

14. An optical system for measuring a diameter, distri- 
bution and so forth of micro gas bubbles and micro 
liquid droplets according to claim 1 3, wherein said 50 
imaging optical system is an anamorphic optical 
system comprising a combination of an axially sym- 
metric objective lens and a cylindrical lens. 

15. An optical system for measuring a diameter, distri- 55 
bution and so forth of micro gas bubbles and micro 
liquid droplets according to claim 13 or 14, wherein 
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FIG. 3 
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FIG. 6 
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FIG. 7 
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FIG. 11 
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